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We consider the radiative decay of atoms scattered by a resonant standing light wave. Scattering 
is shown to suppress the Rabi oscillations and to slow down the atomic radiative decay giving rise 
to a power law behavior of the time-dependent level populations rather than the exponential one. 



INTRODUCTION 

Scattering of atoms by a resonant standing light wave 
is one of the basic phenomena of the atom optics. The 
number of publications on this subject is huge (see, e.g., 
the books ||, [| and the references therein). A very 
special place in this field belongs to a series of works 
on atoms with a wide excited resonance level, the width 
of which is determined mainly by spontaneous radiative 
transitions to nonresonant atomic levels [|[ |[ |6[ ^[ |[ ^) . 
Realizability of such a scheme was demonstrated in the 
experiment || [J. Ar* atoms were prepared initially in 
a metastable state \m) (|1«5, J = 2) = 4s[3/2]2 state, 
correspondingly, in the Paschen and standard spectro- 
scopic notations |)|). In the experiment, the resonant 
light field of the wavelength 801 nm coupled the level E m 
with the excited level E e , which corresponded to the state 
\2ps, J = 2) = 4p[5/2] 2 . 70% of its width was determined 
by transitions to levels different from E m and, finally, to 
the ground level E g (Fig. 1). In principle, in Ar* there is 
another transition, |ls5, J = 2) — > |2p4, J = 2), in which 
98% of the width of the excited level is determined by 
spontaneous decay to levels different from E m ]?J . In any 
case, for both transitions the model of a two-level system 
with a wide excited level decaying predominantly to lev- 
els different from the metastable one works reasonably 
well. This is the model to be considered in this work. 

Both in the experiment || |) and, most often, in the 
theory Q, the investigated regimes of scattering cor- 
responded to the weak-scattering Bragg regime. This 
means that the atomic-beam incidence angle was close 
or equal to the Bragg angle and the resonance coupling 
was not too strong. In terms of the Rabi frequency Q 
and the width of the excited level T, the last assumption 
implies that \Q\ CT. In this case, the decay dynamics of 
an atomic system was shown to obey the usual exponen- 
tial law |t[ |9|. Under some special conditions the effects 
like population trapping were predicted to take place [jrj . 
But, again, the residual atomic population was shown to 



approach its asymptotic non-zero level exponentially . 

In this work we will investigate the dynamics of sponta- 
neous decay in the system under consideration at differ- 
ent conditions. First, we will consider the case of strong 
Rabi coupling, \Q\ ^> T and \Cl\t ^> 1. And, second, we 
will consider the case of normal (or almost normal) in- 
cidence of atoms upon the standing wave. This is the 
diffraction regime of scattering, in which many diffrac- 
tion maxima can arise from the initially well collimated 
atomic beam. 

By investigating the decay dynamics in such a regime 
we find that the total time-dependent populations of both 
metastable and excited levels fall non-exponentially. In 
contrast to standard predictions, the atomic populations 
are characterized by power-law dependencies on the in- 
teraction time t. The effect is not connected with for- 
mation of any kind of grey or dark states (as in Ref. 
0) because asymptotically, at very long time, atomic 
populations tend to zero. But, owing to scattering, the 
radiative decay appears to be slowed down. In addition, 
we find that in atoms scattered by a standing light wave 
the Rabi oscillations of atomic populations appear to be 
strongly suppressed compared to a pure two-level system 
in a resonance field. The physical interpretation of these 
effects is given. 



GENERAL EQUATIONS 

The total wave function \I/ of an atom interacting with 
a light field depends on the atomic center-of-mass posi- 
tion vector r, intra-atomic variables, and time t. The 
wave function obeys the Schrodinger equation 

if = {-^V» + H ot -d.E(r ) t)}* ) (1) 

where Ti = 1, V = <9/<9r, d is the intra-atomic dipole 
moment, and E(r, t) is the electric field strength, for a 
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standing light wave given by 

E(r, t) = 2E cos(c^) cos(kx). (2) 

Here and below x is the center-of-mass coordinate along 
the axis parallel to k; k and Eq are the wave vector 
and field-strength amplitude of one of the two identical 
counter-propagating travelling waves forming a standing 
light wave (Fig. 1), k = uj/c. 



with the matrix Hamiltonian 

it = f-^l ~ ^Po x V x -±ficos(fcc) \ 

{ -^cos(kx) -£-Vl-±p 0x V x -fJ- 

(6) 

Here and below = d/dx, Q = 2d me • E is the Rabi 
frequency, and d me = (ra|d|e) is the dipole matrix ele- 
ment. 



ADIABATIC APPROXIMATION 




Figure 1. A scheme of atom- light scattering and the internal 
structure of atomic levels. 



Explicitly, equations for (p m (x,t) and op e (x,t) equiva- 
lent to (||) have the form (in the case of normal incidence, 
Pox 0): 



ij^(fm(x,t) = - ^Vl(p m (x,t)-^ cos(kx) (p e (x,t) (7) 



d 1 

i — (p e (x,i) = -—V 2 x (p e (x,t) 



d 
and 



With respect to intra-atomic variable, the wave func- 
tion \I/ can be expanded in a series of intra-atomic field- 
free wave functions \i) 



* = y^iQM)|z)exp ( —iEit + ipo 



2m 



(3) 



where po is the unperturbed-atom center-of-mass mo- 
mentum. 

In this work we will consider only the resonance case, 
when the light frequency equals the energy spacing be- 
tween some two discrete nondegenerate atomic levels, E e 
and E m , uo = E e — E m . We assume that E m and E e are, 
respectively, an infinitely narrow metastable and a wide 
excited atomic levels, and the width T of the excited level 
is determined predominantly by its spontaneous decay to 
the ground atomic level (inset of Fig. 1). 

In the resonance case we keep only two terms in the 
expansion (||) with i = m and i = e. Moreover, in the 
rotating wave approximation we retain only one of the 
two terms in the Euler expansion for cosine cos (cot) = 
\ [exp (iout)+exp (—tout)] to drop the fast oscillating terms 
oc exp (±iuut). The arising equations for the metastable- 
and excited-state center-of-mass wave functions (p m (x,t) 
and op e (x,t) can be written in the form of a matrix 
Schrodinger-like equation for the two-component func- 
tion 



$(x,t) 



(ip m (x,t) 1 m 

\ ip e (x,t) J " 



.d<S>(x,t) 
1 di 



H$(a?,t) 



(4) 



(5) 



r n* 

-i-(f e (x, t) - — COs(kx) (p m (x, t). (8) 

In this work we will use the adiabatic approximation in 
which we will drop the kinetic energy operator — V^/2m 
in the Hamiltonian H (||) and, hence, the terms propor- 
tional to — V^/2m on the right-hand side of Eqs. (0) and 
(i to get 



H « H 



ad 





■\£L* cos(kx) 



cos(kx] 

-4r 



(9) 



d Q 

i — ip m (x,t) = -—C08(kx)(f e (x,t), (10) 



and 



i^ip e (x,i) = -i^<p e (x,t) - ^-cos(kx)cp m (x,t). (11) 

For the kinetic energy term (— V x /2m)up rn (x, t) to be 
dropped from Eq. ([?]), it must be smaller than the other 
two terms retained in Eq. (|io|). Hence, qualitatively, the 
applicability criterion for Eq. (|Io| ) has the form 



2m \dt /' 



(12) 



where angular brackets denote averaging over the state 
\I/ . As for Eqs. (||) and (|ll|), they contain additional 
terms ex T, which can be large, and the adiabaticity cri- 
terion for Eq. (|ll| ) is given by 



2m 



<C max < T 



(13) 
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By comparing the conditions of Eqs. ( fL2| ) and (|l3|), we 
see that the first of them is sufficient for the kinetic en- 
ergy operator to be dropped from the Hamiltonian H 
(||) and for the adiabatic approximation to be valid. An 
explicit form of the criterion ( |l2| ) is discussed below. 

If T > (d/dt), the condition (Q can be invalid 
whereas, still, the condition ( |l3| ) can be satisfied. In 
this case the kinetic energy operator in Eq. (Q) has to be 
retained whereas in Eq. (||) it can be dropped, and the 
equations to be solved are (^) and (|ll|). Such a general- 
ization of the adiabatic approximation will be considered 
elsewhere. 

By expressing cp e (x,t) from Eq. (|Io| ) via (p m (x,t), 
tp e (x,t) = —2i(p m /Qcos(kx), and substituting (p e (x,t) 
into Eq. (|ll|), we can reduce the two equations (|iq|) and 
( pj]) to a single second-order equation for (p m (x,t) 



|Q| 2 cos 2 (kx) 



<Pm 0. 



(14) 



As the Hamiltonian 
lutions of the form 



QUASIENERGY SOLUTIONS 

is stationary, Eq. 



<fr(x,t) = exp(— ijt) u 7 (x), 



has so- 



(15) 



where 7 and u 7 are the quasienergies and quasienergy 
wave functions to be found from the eigenvalue equation 



H ad u 7 (x) = ju 7 (x). 



(16) 



In the adiabatic approximation solution of the 
quasienergy problem is very simple. As the operator H a d 
(||) does not contain any derivatives over x, the coordi- 
nate x plays the role of a parameter or a quantum num- 
ber, and the position-dependent eigenvalues of H a d are 
easily found to be given by 



7±0*0 



£2 

T 



\Q\ 2 cos 2 (kx). (17) 



The quasienergy 7+(x) is a direct generalization of the 
Chudesnikov-Yakovlev complex potential j| (at zero de- 
tuning), which follows from Eq. (U7]) in the limit |0| <C T 



, |^| 2 cos 2 (fex) 
7+(ar) « -i — = Vcy(x)\ 



(18) 



In a general case, the complex quasienergies j±(x) ( |l7| ) 
determine both average 'center-of-mass' values Re(7±) 
and widths F± = — 2Im(7±) of the atomic quasienergy 
levels. Explicitly, the widths T±(x) are given by 



T ± (x) 



=F Re 




\ft\ 2 cos 2 (kx) 



(19) 



The broadened quasienergy levels or zones are described 
in Fig. 2. The curves at these picture correspond to 





Figure 2. Quasienergy zones of an atom in a standing light 
wave in the cases 2|fi|/r = l/y/2 (a) and 2^2 (b). 



boundaries of zones determined as Re[7±(#)] + ^T±(x) 
and Re[7±(x)] — ^T±(x). The spacings between the 
boundaries are equal to the widths of the zones T+(x). 
Two different zones are indicated by different shading. 
The pictures (a) and (b) correspond to weak (2\Q\ < T) 
and strong (2|Q| > T) resonance or Rabi coupling of lev- 
els Em and E e . In the case of weak Rabi coupling, one of 
the zones is much narrower than another, and the wide 
zone can be eliminated adiabatically to give rise to the 
description in terms of the potential Vcy(%) ([H)- Math- 
ematically such an elimination of a wide quasienergy 
zone is equivalent to dropping the second-order derivative 
term in Eq. (14), which gives 
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V C Y(x)ip m (x,t). 



(20) 



The condition under which the second-order derivative 
in Eq. ( fL4| ) can be dropped is easily estimated with the 
help of Eq. @: Cp m ~ Vcv^m ~ Vg Y (p m . Hence, 
<p m < T(p m if \V C y\ «Tor |ficos(fca;)| < T. 

In the case of strong Rabi coupling such adiabatic ap- 
proximation and all the resulting equations can be invalid 
at x close to the branching points of the root square in 
Eq. (|l7|). Indeed, at these points the derivative X7 X be- 
comes infinitely large and the kinetic energy terms in 
Eqs. @, (0), (|| ) cannot be dropped. But in the region 
of x close to 7r/2& (which is most important for the given 
below long-time analysis) even in the case |0| ^> T the 
width T+(x) @ is very small and can be approximated 
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by a parabolic dependence on x — w/2k 



I n | 2 cos 2 (kx) \n\ 2 k 2 



r+«-2im(y cy ) 



2k) ■ 
(21) 

As in this region r_ « T ^> T+, we get again a narrow 
quasienergy zone at the background of a wide one, and 
again adiabatic elimination appears to be applicable. 

The role of " nonadiabatic" points, where T « 
2|ficosfor|, will be discussed elsewhere. In principle, the 
arising peculiarities can be observed in experiments with 
scattering of narrow atomic wave packets aimed specifi- 
cally at these points. Such a formulation of the problem 
will be discussed separately too. 



SOLUTION OF THE INITIAL- VALUE PROBLEM 



The found above quasienergies J±(x) (|17|) are suffi- 
cient for solving the initial- value problem. By assuming 
that the interaction is turned on suddenly at t = and 
that <p m (x, 0) = 1 and ip e (x, 0) = 0, we present the time- 
dependent functions (pm ie (x,t) in the form of superposi- 
tions 



1>e (M) = ^4 ±) eWexpH7±Wt), (22) 



where the coefficients A^e(x) are to be found from the 
initial conditions, which yield 



am + a^ = i, 4+) + 4-) = o, 

and the equations following from Eq. 
Hamiltonian (H)) 



2 7± 



Q, cos (kx) 



a(±) 



(23) 
(with the 



(24) 



Eqs. 



and 



and (24) are solved easily to give 



4 ±} = t 



27i 



\/-r 2 + 4|fi| 2 cos 2 (for) 



$7* cos (fcx) 



A /-r 2 + 4|fi| 2 cos 2 (for)' 



(25) 



(26) 



The corresponding time-dependent center-of-mass 
atomic wave functions are given by 



27 T (x) exp(-i^±(x)t) 



MM) v /_ r2 +4 |n| 2 cos 2 (^) 



and 



^* cos(/cx) exp(— ij±(x)t) 
v/-r 2 + 4|^| 2 cos 2 (te) 



The squared absolute values of the functions (p mi e(x,t) 
determine the probability densities to find an atom at a 
time t in a vicinity of a point x at the levels E m and E e 

- \(Pm ie (x,t)\ . (29) 

CtX 7T 

Integrated over x from zero to 7r/fc, the probability densi- 
ties dW m , e (x, t)/dx give the time-dependent total proba- 
bilities of scattering at a single period of a standing light 
wave 



"V.r'm = / dx 

Jo 



V/ \ dW m , e (x,t) 



dx 



7T 



dx | <p mi e(x,t) | 5 



(30) 



LONG-TIME ASYMPTOTIC LIMIT 



The long-time asymptotic limit corresponds to the case 
Ft ^> 1. In this case the main contribution to the integral 
over x in Eq. ( |30| ) is given by the most slowly decaying 
terms. In the case of a strong Rabi coupling \Q\ ^> T 
such slowly decaying terms correspond to the quasienergy 
7 + (x) and to the region of x close to n/2k. By assum- 
ing that in this region the product \ft\ 2 cos 2 (kx) is small 
compared to T 2 /4, we can reduce Eqs. (^?j) and (^8|) to 
the form 



\(Pm(x,t)f 



' exp • 



cos 2 (kx) 



(31) 



and 



\ip e (x,t)\' 1 



cos 2 (kx) exp • 



\n\H 



■ cos 



(kx) 



(32) 

These equations follows also from Eq. (|14|), which is 
solved easily and which is valid under the condition 
\Qcos(kx)\ < 2T. The position-dependent metastable- 
state probability density was found earlier (Eq. (19) of 
Ref. [|| in which the spontaneous decay rate 7 should be 
substituted by |^| 2 /2T). 

At x close to 7r/2/c and with cos (kx) approximated by 
the parabolic function (^l|), Eqs. (31) and ( |32| ) take the 
form 



\<Pm(x,t)\' 



^{Jip( kI -iy-} (33, 



(27) and 



\(f e (x,t)\ 2 



^(--I) 2 -{-^(--|) 2 }- 

(34) 

(28) The functions \(p m (x,t)\ 2 (33) and \(p e (x,t)\ 2 (Q) are 



plotted in Fig. 3. The width Ax of the interval where 
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Figure 3. Probability densities to find an atom at the levels 
E m and E e around a position x; \Q\/F = 3, Ft = 2. 



\(p m (x,t)\ 2 and \(p e (x,t)\ 2 are not small is given by 



Ax 



1 



k\n\ 



The condition Ax Cl/fc has the form 

\n\H 1 



(35) 



(36) 



This is the condition under which only relatively narrow 
regions close to nodes of a standing light wave give not 
small contributions to the total probabilities of scattering 

w m ,S) ©• 

The parameter Ax ( |35| ) is one of the key parame- 
ters of the problem under consideration. In particu- 
lar, the condition Ax <C 1/k (or the condition given 
by Eq. (|3rj)) justifies the approximation of cos 2 kx by 
a parabolic function in Eqs. (^3|), and (|34|). As for 

the used above approximation |^| 2 cos 2 (kx) <C T 2 /4 for 
\x — 7r/2k\ < Ax, its validity is related to both assump- 
tions, Ax <C 1/fcandrt ^> 1. Indeed, at \x— ir/2k\ ~ Ax, 
the "effective" position-dependent Rabi frequency 



Q e ff(x) = Flcos(kx) 



is estimated as 



effl 



\Q\kAx 



(37) 



(38) 



The ratio \Q eff \/T ~ 1/VTt is small if Ft > 1 whereas 
in the opposite case, Ft < 1, |f2 e //| > T (and, of course, 
O e //^ ^> 1). Hence, we expect that transitions between 
the metastable and excited levels have significantly dif- 
ferent form at long and short times, Ft ^> 1 and Ft < 1. 
In the first of these two cases (long-time asymptotic) the 
transitions E m — > have a form of irreversible transi- 
tions to the quasicontinuum of the wide excited level E ei 
whereas in the second case (short-time limit) they take a 
form of multiple Rabi oscillations. In the following sec- 
tion we will see how the dynamics of excitation is affected 
by a mixture of these two types of transitions. 



The parameter Ax ( |35| ) determines the relation be- 
tween the heights of the curves | ip m (x, t) | 2 and | cp e (x, t) | 2 : 

if \Vm\max = ^ 



2 _ 1_L^ 2 At 2 

e\max ~ p2 



Ft 



«1, 



(39) 



if Ft > 1. 

As Ax ( |35|) is the width of the " most important" region 
of x, where T + (|2|) is small and the corresponding part 
of atoms decays slowly, we can estimate now a rigidity 
of the assumption about the normal incidence of atoms 
upon a standing wave. If po x 7^ 0, the atoms move ho- 
mogeneously along the x-axis. The arising displacement 
during the interaction time is po x t/m, and it must be 
not larger than Ax to keep atoms decaying slowly, which 
gives 



fc| ^ l ^MVT <<1 ' 



(40) 



where vo x = Pox/m, and the last inequality follows from 
Eq. ©. 

At last, the definition of the characteristic "important" 
interval Ax ([35]) can be used to evaluate the validity 
criterion of the adiabatic approximation. In this ap- 
proximation the characteristic value of the kinetic energy 
— V 2 /2m is assumed to be small compared to the charac- 
teristic value of d/dt. The latter is estimated as (at Ft ^> 
1 and |x - ir/2k\ ~ Ax): d/dt ~ T+ ~ |ft| 2 /c 2 Ax 2 /r ~ 
1/t. In accordance with the uncertainty principle, we 
put Va; ~ 1/Ax, which gives V 2 /2m ~ uu r \Q\ 2 t/F, where 
uo r = k 2 /2m is the recoil frequency. With the help of 
these estimates the applicability criterion of the adiabatic 
approximation V 2 /2m <Cl/t can be reduced to the form 



(41) 



As, typically, oj r ~ 10 3 r, this condition can be fulfilled 
both at \tt\t < 1 and \Q\t > 1. 

By returning to the analysis of the time behavior of the 
scattering probabilities, note first that the decay of small 
portions of the atomic wave function localized near any 
given x has an exponential character, though with the de- 
cay rate depending on x (Eqs. ( |33| ) and (34)). But, as the 
sum of exponents is not identical to any other exponential 
function, it's not surprising that the total probabilities of 
scattering ( ^0|) decay non-exponent ially. In other words, 
the total probabilities of scattering W^ e \t) are deter- 
mined by areas under the curves |^ m , e (x, t)\ 2 at Fig. 3. 
With a growing time t these areas shrink, but the laws 
of their decreasing are not exponential. 

In principle, in the asymptotic limit Ft ^> 1 the inte- 
grals in ( pQ| ) are easily calculated with cp m ^ e (x,t) substi- 
tuted from Eqs. ( |33| ) and ( |34| ) and the limits of integra- 
tion extended to ^00. But it's very interesting and in- 
structive to use Eq. (^T|) before calculations to reduce the 
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integrals over x to integrals over the narrow-zone width 



r+: 



Vr r°° dr + 



-r+t 



7T|fi| Jo Jr^ 



pl/2 



(42) 



and 



oo 

/ 



r + dr + e 



1 



2|^iv^rt 3 / 2 



(43) 

So, indeed, the long-time behavior of the total probabil- 
ities to find an atom after scattering at the metastable 
and excited levels is determined by power-law rather than 
exponential dependencies on the interaction time t. 

It should be noted, that, in principle, the non- 
exponential decay characterized by Eqs. (f|2|), (|i^ ) can 
occur also in the case of weak Rabi coupling, |0| < T, if 
only the conditions j3^), (^5|) , and jil] ) are fulfilled. The 
first of these conditions (|3q ) at |0| < F can be fulfilled 
only if the interaction time t is very large, t ^> T/|^| 2 ^> 
1/|^| ^> 1/r. Such a long time can make the restriction 
of the transverse velocity vq x ( |4C| ) too severe to be easily 
satisfied. For this reason, the case of strong Rabi cou- 
pling looks much more favorable than the case of weak 
coupling for observation of the effects described in this 
and the following sections. Compatibility of the condi- 
tions @ and @ requires the Rabi frequency to be not 
too small, \Vt\^> vWr. 



PARTIAL PROBABILITIES OF SCATTERING 
INTO DIFFRACTION BEAMS 



These equations look like equations for two coupled an- 
harmonic oscillators in a resonance field with the reso- 
nance detuning 5 — kpo x /m and the anharmonicity pa- 
rameter coinciding with the recoil frequency uo r . The 
above-discussed adiabatic approximation corresponds to 
ignoring the anharmonicity terms n 2 u r in Eqs. (|46|). 
With the additional assumption about the normal inci- 
dence, pox = 0, Eqs. (pS) take the form 



• (m) 



(*) = -°(«£i+«#i). 



ia 



(e) 



(*) = 



iF 



(e) 



ft* 

T 



(a^i + 'ffi)- (47) 



The transition from Eqs. ( |46| ) to ( |47| ) can be consid- 
ered as the Raman-Nath approximation for the two- 
dimensional (2D) system. It should be noted, however 
that both Eqs. ( f47| ) and their solutions presented below 
differ significantly from and are much more complicated 
than the standard ID Raman-Nath equation and its so- 
lution [@. 

The simplest way of finding a^ 1 ^ (t) obeying Eqs. (47) 
is related to the calculation of the Fourier transforms (44) 
of the earlier found functions (p m (x,t) ( |2^ ) and (p e (x,t) 
(f27|). Not dwelling upon the details of calculations, let 
us present here the arising results: 



J2n 



\Q\t 



Ft 

T 



/* 

Jo 



dzJ2n 



\n\t 



To investigate in more details the time evolution of 
scattering, let us consider the Fourier transforms of the 
atomic center-of-mass wave functions (p mi e(x,i) 



(m,e) 



_ k r 

Jo 



dx ip m ,e ( x , t) exp (—inx) . (44) 



The functions a^' e \i) and their squared absolute values 



and 



■A 



(48) 



(m,e) 



(*) 



(45) 



are the probability amplitudes and partial probabilities 
to find an atom at the levels E m or E e in the n-th diffrac- 
tion beam with the momentum po + nk, which makes an 
angle n « nk/po with po, n = 0,±1,±2, .... Equations 
for a^™' e \i), equivalent to Eqs. (||), (^), are given by 

ia$T\t) = («V + nS) - £ («£>! + 0&) , 



J: 



2n+2 



whereas ^^(t) = and W£'(t) = 0. In Eqs. (|4j), 



^2n 



00/ 



\n\t 



(49) 



J2 n are the Bessel functions and Iq and I\ are the 
modified Bessel functions [11]. It should be emphasized 
that in derivation of Eqs. (48) and ([49|) we did not make 
any assumptions about a value of the parameter Ft and, 
hence, these equations are valid both at short interaction 
times, Ft < 1, and in the long-time asymptotic limit, 



Ft > 1. The functions W$™\t) 



8\) and^^t) 



) 



(46) are plotted in Fig. 4a and b. In contrast to the 
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b 

Figure 4. The time- dependent partial probabilities Wn m \t) 
(a) and W^ e) (£) (b); |ft|/r = 5. 

well-known Rabi oscillations in a pure two-level system 
in a resonance field, oscillations of (t) are aperiodic 
and both positions of peaks and zeros of Wc^^ (t) de- 
pend on n. 

Alternatively to (|30|), the total probabilities of finding 
atoms after scattering at the levels E m and E e can be 
defined as sums over the diffraction beams 

w4rV) = Ei a ^ e) wi 2 - ( 5 °) 

n 

The dependencies W$%\t) and W^]{t) © calculated 
numerically by summing W^\t) © and W^ +1 (t) © 
are shown in Figs. 5a and b. In contrast to W^\i) and 

W^n+iit), oscillations of W^\t) and W^}(t) are peri- 
odic, and their period coincides with that of the Rabi os- 
cillations in a pure two-level system driven by a resonant 
field with the field-strength amplitude 2Eo (the dashed 
curves in Fig. 5) . However, as it's seen well from Fig. 
5, the quasi- Rabi oscillations of the functions W^\t) 



and Wfol (t) for scattered atoms are strongly suppressed 
compared to those of a two-level system. The effect of 
suppression of Rabi oscillation is explained mainly by a 
kind of inhomogeneous broadening. Oscillations of par- 
tial probabilities w£\t) © and W^ +1 (t) © are well 
pronounced and their amplitudes are large enough (Fig. 
4). However, as "periods" of these oscillations are differ- 
ent for different n, summation over n smoothes over the 
oscillations and decreases their amplitudes. In insets of 




0.5 1 1.5 2 2.5 3 4 



Figure 5. Total populations W^\t) (a) and W^ t {t) (b) at 
|0|/T = 5 (solid lines) and the pure two- level-system popula- 
tion (dashed curve) 

Fig. 5 it's seen clearly that the degree of suppression 
the Rabi oscillations increases with a growing value of 
Tt: in the case Tt ^> 1 oscillations are almost completely 
smoothed away. This effect is related to the discussed 
above change in the character of transitions between the 
levels E m and E e : Rabi oscillations in the case Tt < 1 
and "discrete level - quasicontinuum" transitions in the 
case Tt ^> 1. 

Another effect which is clearly seen in the insets of 
Figs. 5a and b is the scattering-induced slowing down of 
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the radiative decay. The long-time behavior of the prob- 
abilities W^ e \t) is described rather well by asymptotic 



formulas of Eqs. (|42|) and (|43|) . The difference with ex- 
ponential decay in a pure two-level system is rather well 
pronounced. 



CONCLUSION 

So, the main predictions of the carried out considera- 
tion are (i) suppression of the Rabi oscillations in the case 
of atom scattering (compared to a pure two-level system) 
and (ii) slowing down the radiative decay and forma- 
tion of nonexponential (power-law) tails in the dependen- 
cies W^ e \t). Qualitatively, our interpretation of these 
effects consists of the following. The Rabi oscillations 
are suppressed owing to inhomogeneous-broadening-like 
effects when the partial probabilities of scattering are 
summed over the diffractions beams. The power-law 
dependencies and slowing down of the radiative decay 
arise because of the position-dependent modulation of 
the field-strength amplitude in a standing light wave 
and, hence, the position-dependent modulation of the 
decay rate of slowly decaying quasienergy atomic levels 
T+(x). As the result, of this modulation the manifold of 
the arising quasienergy levels is characterized by contin- 
uously varying with x and approaching zero at x = ir/2k 
width. Populations at these quasienergy levels decrease 
exponentially but with different, x-dependent rates, and 
their superposition gives rise to the non-exponential de- 
cay laws. 

We assume that observation of these effects can be 
made in the framework of an experiment similar to |], 
though with some modifications, to provide the con- 
ditions for the strong Rabi coupling and the normal- 
incidence diffraction regime of scattering. To compare di- 
rectly the time evolution of atomic populations in atoms 
scattered by a standing light wave and in a pure two-level 
system, one can make two series of similar measurements: 
in a standing light wave and in a single travelling wave 
of a doubled field strength amplitude. In the last case 
the standing-wave scattering effects disappear and the 
large-amplitude Rabi oscillations and the usual exponen- 
tial decay have to be observed. 

Finally, the above-described scattering-induced sup- 
pression of the radiative decay in the case of strong Rabi 
coupling reminds the effect of interference stabilization in 
Rydberg atoms in a strong light field [[l2|, |l3| ]. In both 
cases the effects of slowing down the decay processes are 
related to formation of narrow quasienergy levels and in- 
terference of transitions from different levels to the com- 
mon continuum. We find this analogy important because 
it establishes links between different regions of physical 
phenomena and demonstrates a rather general character 
and fruitfulness of the idea of interference stimulated by 
sufficiently strong interactions. 
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